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1 Introduction 



Owing to the observation of quasi two dimensional antiferromagnetism in the undoped 
phase of high T c superconductors 0, the theory of two dimensional quantum antiferro- 
magnetism has received much attention recently. The continuing discussion focusses on 
the question of whether or to which extent the classical picture of a Neel type ordered 
ground state remains valid in the presence of quantum fluctuations which are known to 
be particularly important in low dimensional systems. For the spin-1/2 isotropic Heisen- 
berg antiferromagnet (HAF) with nearest neighbour interaction on a square lattice this 
question has been settled: the result of numerous quantum Monte Carlo (MC) simulations 
leave no doubt about the existence of Neel order in the ground state of the square 
lattice HAF. Moreover, by refined MC techniques it has been possible to determine the 
staggered magnetisation M. of this state with high accuracy. Quantum fluctuations are 
found to reduce M. to about 60% of its classical value, but surprisingly the relative differ- 
ence between the most accurate numerical value M. = 0.3074(4) and the value obtained 
in first order spin-wave (SW) approximation Adx.sw — 0.30340 is only about 1%. Long 
range order is opposed not only by quantum fluctuations but also by frustration. For a 
spin-1/2 Heisenberg model, frustration is defined in the sense introduced by Toulouse Q 
for the z-z part of the coupling and, with respect to the x-x and y-y parts of the inter- 
action, by the impossibility of using the Marshall Ansatz |5[ for the ground state wave 
function ||. The HAF on a triangular lattice (HAFT) which will be studied in this work 
is the prototypical example of a frustrated antiferromagnet. In early work, Anderson and 
Fazekas J/J suggested that the combined effect of quantum fluctuations and of frustration 
would favour a disordered overall singlet state over an ordered Neel type ground state. 
The state proposed by Anderson and Fazekas consists of products of singlet pair states 
and is commonly referred to as a resonating valence bond (RVB) state. It has received 
much attention in the theory of high T c superconductivity. Numerous subsequent investi- 
gations have led to contradictory views of the ground state of the HAFT. The spin-wave 
analysis || [| predicts that quantum fluctuations are insufficient to suppress the classical 
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Neel order. Certain variational calculations support this conclusion |TIJ, others contradict 
it. In fact, an exotic RVB-type state which breaks time reversal invariance and parity 



has also been obtained by a variational method Several attempts have been made 

to determine the nature of the ground state of the HAFT by extrapolating from results 
obtained by the exact numerical diagonalisation of the Hamiltonian of small periodic cells 



12] , y|. In all of these studies the results of the extrapolations have been interpreted as 



being indicatory of a disordered ground state or of a state that is close to the critical point 
of losing long range order. The analysis of a high-order perturbation expansion has lent 



support to this view |[14|| . However, by extending the numerical studies of finite cells to 
include not only the ground state energy but the entire spectrum of the HAFT and by 
an appropriate interpretation of their results Bernu et. al. \IE\ have very recently given 
convincing arguments in favour of a Neel ordered ground state of the HAFT. In view of 
the experience with spin-1/2 square lattice HAF, it might thus be hoped that the SW 
approximation also yields accurate results for the physical quantities of the HAFT. It is 
one of the purposes of this study to show that the results obtained by the exact numerical 
diagonalisation of small cells of the HAFT do indeed point towards the validity of the SW 
approximation. 

With regard to this purpose it will prove useful to consider the HAFT with nearest 
(nn) and next nearest neightbour (nnn) interaction, 

N N 

h = j2 s Sj + a J2 s i s J ( L1 ) 

nn nnn 

where Sj denotes a spin-1/2 operator at the lattice site R,. (We set the nn coupling equal 
to unity without loss of generality). As is seen in Fig. 1, the nnn interaction couples spin 
pairs within each of the three sublattices. Therefore, for sufficiently strong ferromagnetic 
nnn coupling, —a ^> 1, the spins of each sublattice will order ferromagnetically, forming 
a macroscopic spin. In this situation, a SW expansion will certainly be applicable. Thus, 



the parameter region —a ^> 1 of the model provides a firm basis for a comparsion of 
the SW approximation and of the exact numerical result. 

The model ( |1 . 1| ) is also of considerable interest in its own right. In the classical approx- 



imation one finds that its ground state changes discontinuously from the 120° three sublat- 
tice state for a < a c i = 1/8 to a four sublattice state for a c i < a < 1 |16| , Fig. 2. The four- 



sublattice state is degenerate with respect to the angle 9, i.e. its energy is E = —2(1 + a) 



independent of 9. As has been shown in Ref. [0 quantum fluctuations select the state 
with 9 = from the continuum of degenerate four sublattice states. We shall refer to this 
state, which is a two sublattice structure, as the collinear state. At a = 1 the ground state 
structure changes again to become incommensurate for a > 1. In the limit a — > oo, the 
three sublattices decouple so that a 120° structure exists on each of them. The manifold 
of classical ground states of (|1 . 1|) that evolves as a increases from — oo to oo is concisely 
described by 

Sj = ei cos QRj + e 2 sin QRj, (1.2) 

where e± and e 2 are an arbitary pair of orthogonal unit vectors and where Q = Q(a) traces 
out the path in the Brillouin zone shown in Fig. 3. 

The question of how the transition between the 120° and the collinear state of the HAFT 
with nn and nnn couplings is affected by quantum fluctuations has attracted considerable 
interest [JTBL |T7], [TS], [BJ. If both types of long range order survive at the quantum level, 
the transition between them may be of first order, i.e. the order parameters of both states 
may remain finite at the transition point so that the symmetry breaking pattern changes 
discontinuously. Alternatively, one or both types of long range order could be suppressed in 
the vicinity of the classical transition point, and an intermediate ordered or disordered state 
could appear at the quantum level. Disregarding the first possibility, Baskaran conjectured 
that a time reversal and parity breaking chiral spin liquid state should evolve as the ground 
state of the model (|1.1|) when a exceeds a critical value [|17]]. In a previous study Jl8| 



we 



could not confirm this conjecture. More recently, Chubukov and Jolicceur []1{J investigated 
the collinear structure by a self-consistent SW approximation. These authors find that 
the stability region of the collinear structure extends to values of a < a c i and argue that 
the transition between the 120° and the collinear structure should be of first order. It 
is the second purpose of this paper to provide and discuss further results concerning this 



transition. 



In section 2 we work out the second order SW expansion for physical properties of the 
HAFT pertaining to the 120° structure, and we show how finite size results are obtained 
within the SW approximation. Section 3 contains our numerical results for physical proper- 
ties that characterise the ordered states of our model ( p..l|) : the sublattice magnetisations, 
the helicity and the chirality. The emphasis is on a direct comparison of the exact numeri- 
cal results with the predictions of the SW-theory for finite systems and on the discussion of 
the behaviour of the order parameters in the transition region between the 120° structure 
and the collinear structure. In section 4 we summarise our conclusions. Technical details 
are deferred to three appendices. 

2 Spin-wave approximation 

To order 1/S, the spin- wave approximation for the model ( |1 . 1|) has been discussed by 
Jolicceur et al || and a calculation of the 1/S 2 corrections to the spin- wave spectrum of 
the collinear state has been performed by Chubukov and Jolicceur fllDfl . Furthermore, for 
a = the 1/S 2 correction to the sublattice magnetisation has been obtained by Miyake 
P|. Below we present the results of an extension of Miyake's calculation to the range 
— oo > a > 1/8 where the 120° state is supposedly stable. With the objective of comparing 
the SW approximation with our numerical results we also include the first and second order 
results for the helicity of the 120° structure, and we discuss the finite size corrections to 
the sublattice magnetisation and the helicity. 

The starting point of the SW expansion is the Hamiltonian 



H 



nn nnn 



cos% (S?S? + S*S*) + sin% (SfS? - S*S*) + SfS] 



(2.1) 



with 6 a 




which is obtained from ( |1.1| ) by rotating the local z-axis to the direction of the classical 
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sublattice magnetisation: 

S? = caa9 i Sf + fnn9 i S;, = Sf, S? = cos ft Sf - sin ft Sf. (2.2) 

For the 120° state ft = 0,27r/3, — 27r/3, when i is a site of the sublattice A, B or C while 
ft = 0, 7r on the sites of the two sublattice A or B of the collinear state. B represents 
a staggered magnetic field. From (|2.1|) the SW expansion is obtained by replacing the 
spin-operators Sf by Holstein-Primakoff boson operators and by expanding with respect 
to 1/S. In the case of the 120° state the result is 

H = S 2 (H + H 1 + H 3/2 + H 2 ), (2.3) 



with 



H = -^(l+2a)N-^N, (2. 4a) 

Hi = ~ Yl {-\( a t a J + h - c -) + ^( a t a t + h - c -) ~ \( n i + n j)} 

+^T, n i + ^J2{( a t a i + h - c -)-( n i + n j)}^ ( 2 - 4b ) 



H3/2 



: -J-^J2 Sind ij( n i a 3+ h - C -)> {2. Ac) 

V & nn 

H 2 = -—Yl-{nin j ) + -(n i a i af + h.c.)--(n i a i a j + h.c.)\ 

o nn 12 J 

+ C2 X, \ nin 3 ~ ( n i a i a t + h - C -)\ ■ ( 2 - 4d ) 



Diagonalisation of Hi yields the SW frequency (we adopt the notation of Miyake || wher- 
ever possible): 

u(k,B) = 3Su(k,B) 



3,S'\/ ( L + ^ - 2a(k) - 7 (k)) (l + ^ - 2a(k) + 2 7 (k) ) . (2.0) 



where 

a 



a(k) = a - ^ [cos(V3k y ) + cos(^ + ^) + cos(^ - ^)) , (2. 6a) 
7 (k) = i(cos(M + cos(| + ^)+cos(|-^)). (2.6b) 



For a < 1/8, u;(k, 0) = uj(k) vanishes linearly in the centre of the Brillouin zone (BZ) and at 
opposite corners k = ±Qi = ±(47r/3, 0). The zero modes are the three Goldstone modes of 
the symmetry broken state f2~0fl . The two modes at ±Qi, which are degenerate, correspond 
to an infinitesimal tilt of adjacent plaquettes of three spins around two orthogonal axes in 
the plane defined by the ideally ordered classical state. The third mode at k = whose 
SW velocity is larger than that of the previous ones by a factor of v2 corresponds to an 
infinitesimal twist of adjacent plaquettes around an axis that is perpendicular to the plane 
defined by the 120° state. 

As a approaches the value 1/8 from below, the SW frequency softens at Q2 , an d 
Q2 1 Fig- 3, yielding new linear excitation branches at these points. Within this first order 
approximation, these new modes indicate the onset of collinear order at a = 1/8. They 
are, however, unphysical since they are not Goldstone modes of the 120° state. 

In first order in 1/S the sublattice magnetisation is given by the expression |jj 



M l . sw = S jl + -L 1 £(2 - A 1 (k) - A 2 (k))| , 



(2.7) 



where 



A 1 (k) 



l-2«(k) + 2 T (k) Mk) = Mkr . 



\ l-2a(k)- 7 (k) 

To obtain the second order correction M2.SW we follow the method of Miyake and calculate 
the 1/S 2 correction to the ground state energy 

P (2)/ m _ q2, tt v , n2V l(W}l#3/2|{0})|' 

E (a, B) - S (H 2 ) + S £ E(m _ E({nk] y (2-8) 
where {n^} is an eigenstate of Hi and E^n^}) is the corresponding energy 

«) = E(«k+^(k). (2.9) 

The explicit expressions are given in the App. A. As in the case a = considered by 

(2) 1 — 
Miyake, both contributions to Eq (a, B) contain singular parts proportional to y B which, 

however, can be shown to cancel in the sum (12. 81). This feature has been utilised to control 



6 



(2) 

the precision of the numerical computation of the integrals contributing to E$ (a, B) and 

(2,) 

to extract the second order correction M.2.SW to the magnetisation from Eq (a, B) in the 
small B limit, 

_ 1 EP(a,B)-El?\a,0) 
M 2 . sw - hm . (2.10) 

In Fig. 4a we show plots of JM.\.sw an d of the sum M-i.sw + M-2.SW- A plot of M.2.SW is 
included in Fig. 8. Owing to the unphysical zeros which develop in uj(k) at Q 2 *\ Q2 > Q2 
as a approaches a c i the second order contribution M.2.SW diverges at a = a c i. However, 
•M-2.SW enhances the sublattice magnetisation also for a<0.05 where the influence of the 
divergence is still small. Thus, it appears that the stability region of the 120° structure 
extends beyond a c i and hence overlaps with the stability region of the collinear structure. 
This corroborates the argument of Chubukov and Jolicceur according to which the 
transition between the two structures should be of first order. 

In the classical 120° state an Ising-like variable can be attached to each triangular 
plaquette of the lattice [ETJ . It takes the values ±1 depending on whether the spins on the 
three sites of the plaquette turn clockwise or counterclockwise. For the quantum system 
the operator associated with this variable - we shall call it helicity - has been defined by 

X A = (S< x + S, x S fc + S fc x S,,), (2.11) 

where i,j and k denote the sites of an elementary triangular plaquette in a clockwise sense. 
The normalisation has been chosen such that each component of x A has eigenvalues +1,0 



and —1. Long range helicity order is measured by the expectation value [22 



^ = ^E<Xa>, (2-12) 



N A 



where the sum extends over all upward pointing triangles of the lattice. Long range helicity 
order can exist in the absence of long range 120° order. In that case it is the signature 
of a spin nematic [p3 j. It must necessarily exist for a state that possesses long range 120° 
order. In the local reference frame (|2.2p of the 120° state, Sj x Sj takes the form 

(Si x Sj) ■ e y = sin % (SfSf + S?Sj) + cos % (S?S? - S*S*). (2. 13) 



From this expression one immediately obtains the order parameter x v i n first order SW 
approximation 

Xls w {a) = 35 2 |l + ^lx:(2-^i(k)-A 2 (k)-5 1 (k))| 

= 35 2 {l + l{M 1 . sw -S) + ^{StSf) l . sw y (2.14) 



where 

S 1 (k)= 7 (k)^ 1 (k). (2.15) 

A plot of Xi.sw{ a ) is contained in Fig. 4b. As is described in App. A, the second order 
correction X2.sw can be obtained by the same technique as has been applied in calculating 
M 2.sw- A plot of X2.sw{®) is included in Fig. 10. For a< — 1, the quantum fluctuations 
are seen to enhance x y slightly over the classical value. Quite generally, the quantum 
fluctuations are seen to be less effective in reducing the helicity x y than in reducing the 
magnetisation M.. This is not unexpected, since the long wavelength modes that twist 
the spin configuration around the axis perpendicular to the plane of the 120° structure are 
ineffective in reducing x v ■ 

As we have explained above we need to evaluate the SW expressions ( |2.7|) and (|2.14 ) for 



finite periodic systems. This requires modifications, since in ( |2.7| ) and ( |2.14| ) the members 



of the k sums that correspond to the zero modes k = and k = ±Qi are infinite. In his 



seminal work on quantum antiferromagnetism, Anderson |24j has pointed out that these 
infinities are unphysical, yet unavoidable consequences of the basic assumptions underlying 
the S W theory. In previous work on finite size effects [25| it has been shown that the leading 



size correction to Aii.swi®, N) is obtained by omitting the infinite contributions Ai(k = 0) 
and ^b(k = ±Qi) from the sum in (|2.7f ). As we show in App. B this omission yields an 
unphysical result in the limit of strong ferromagnetic nnn coupling, —a/N 3> 1, 

J^ M ^^ = S { 1 + M} >S - (2 - 16) 
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At the same time we find that the expression 

M l . sw (a,N) = S{l + ^± J2j2-A 1 (k)-A 2 (k))\, (2.17) 

in which not only the infinite parts of the zero mode contributions but also the finite parts 
have been neglected, yields the correct result 

lim Mx.sw{oc,N) = S. (2.18) 

— a/N— *oo 



Furthermore, for finite a the leading 1/yN size correction remains unaffected by the 
neglect of the finite parts of the zero mode contributions. Thus, ( |2.17| ) qualifies as a finite 
size approximant for M.\.sw m the entire region — oo<«<1/8 where the 120° structure is 



classically stable. Of course, for a specified finite periodic system of size N the sum ( 2.17 ) 
has to be extended over that discrete set of inequivalent k-vectors which represents the 
system in the k space. Plots of Mi.swi^, N) obtained in this way and of the asymptotics 

1 215 

M 

v iV 



are shown in Fig. 5. Obviously, for iV<10 2 the subleading N~ l corrections preclude the 
use of the asymptotics ( |2.19| ) in an extrapolation to the bulk limit. Fig. 5 also contains a 
plot of 

1215 3 

AW(0, N) = Mlsw^, oo) + —=■ - — (2.20) 

a/ Jy Zl\ 

in which the last term is the subleading correction that obtains in the limit —a/N 3> 1, 
App. B. Evidently, ( p.20| ) provides a good approximation for M.\.sw{^i N). The remaining 
fluctuations in the data points reflect the dependence of A4i.sw(0, N) on the shape of the 
finite cells, i.e. on the specific set of k points to be used in (|2.17|) . 

By analogy with ( |2.17| ) we define the finite size approximants for the helicity by 



xlswi^N) = 3S 2 {l + ^± E (2-A 1 (k)-A 2 (k) + 5 1 (k)) 



• (2.21; 



In principle, finite size values for the second order quantities M.2.SW an d X2.sw can also be 
obtained by omitting the contributions of the zero modes. However, in practice the results 
for Eq \a, B, N) turn out to depend irregularly on \[B so that the extrapolation (|2.10|) is 
not feasible. 



3 Numerical results 



Presently, exact results for the HAFT are available for system sizes up to N = 36 JT3], |15] . 
Thus, in view of the results of the last section there is little use in attempting to extrapolate 
from such system sizes to the bulk properties of the HAFT. Clearly, however, one expects 
to find agreement between the exact results and finite size SW results of the previous 
section, at least in the region —a ^> 1 where the SW approximation has a firm basis. 

In the following subsections we shall discuss the order parameters A4 and M of the 120° 
structure and of the collinear structure (subsections 3.1a and 3.1b). The helicity x an d the 
chirality P will be studied in subsections 3.2 and 3.3. The numerical calculations include 
cells of sizes iV = 12, 18, 21, 24, 27 and 30 with periodic boundary conditions, Figs. 6a-f. 

The translationally invariant N = 21 and 27 cells are fully symmetric under the C§ v 
symmetry group of the Hamiltonian but they do not accommodate the 2x1 cell of the 
order parameter of the collinear state. Since we are interested in obtaining both the order 
parameter of the 120° state, whose unit cell is of size \/3 x v3, and that of the collinear 
phase we include the N = 18, 24 and 30 cells, which accommodate both unit cells but are 



not symmetric under all C% v operations. We diagonalise the Hamiltonian (|1.1 ) for these 
finite systems by using the Lanczos algorithm. For a < 1/8, the ground state of cells 
with an even number of sites are translationally invariant singlets, S to t — 0, Q = Qo = 0, 
while the ground state of cells with an odd number sites are doublets, Stot — 1/2, with 
wavevector Qi = (47r/3,0). At some value a > 1/8 the ground state of the iV = 18 and 
the N = 30 cells cross over to the wavevector Q 2 , which is the Bragg vector of the collinear 
state. As has been discussed in [113 for the case N = 21, when a exceeds a critical value 



&q U {N) > 1/8 (a gu (21) = 0.253, a qu (27) = 0.130) the ground state of cells with odd N cross 
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over to a wavevector Q 7^ Q2, since the periodicity of the collinear state is incompatible 
with the periodicity of these cells. We notice that the same type of behaviour already 
occurs at the classical level. For the N = 21 and the N = 27 cell the classical ground state 
(|0|) jumps from Q = Q : to a Q ^ Q 2 at a cl (21) = 0.173 and a cl {27) = 0.129. 



3.1 Sublattice magnetisation 

3.1a The 120° order 

As a measure for the sublattice magnetisation Ai of the 120° structure we calculate the 
structure function 

N 

MQi)=E elQl(ri ~ rj) ( S * S ;>o with Qi = (4tt/3,0) (3.1) 

from our numerically determined ground states. To find the connection between SV(Qi) 
and the finite size approximants A4(a, N) for the sublattice magnetisation we consider the 
limit of large ferromagnetic nnn coupling, —a 3> 1. In this limit the spins of each of the 
sublattices A, B and C will align ferromagnetically forming macroscopic spins of magnitude 
Sa = Sb — Sq — N/6. These macroscopic spins may fluctuate against each other so 
that the total spin S to t = |Sa + Sb + Sc| agrees with the values found in the numerical 
calculation, i.e. S to t(N) = for even iV and S to t(N) = 1/2 for odd N. A simple calculation 
which is detailed in App. C yields 

lim = ( * + * ' JVC ™' (3.2) 



Since 



it follows from (|3.2j ) that in the limit a — > — 00 



1 2 (^r 11 - iiv) , A^even, 

M\a,N) = \ ) N q iN > ' (3.3) 
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Here, the factor of 2 reflects the fact that SW(Qi) is the structure factor of only one of 
the two 120° structures with wave vectors ±Qi which are contained in the ground state 
with equal weight. The above considerations suggest that ( |3.3|) still holds approximately 
for — a>l. The same argument also suggests that for — a>l the numerically exact results 
M(a,N) should be well approximated by the finite size SW results Msw( a i N). Plots of 
M\ sw (a, N), i.e. of the first order SW result (|2.17|) , and of M 2 (a, N) as functions of a 
are shown in Figs. 7a-f. The difference 

AM(a,N) = M(a,N)-Mi.sw(a,N) (3.4) 

is displayed in Fig. 8 together with the second order SW correction J\A 2 .sw- Obviously, 
the asymmetric N = 18 cell yields exceptionally large values for AM. Apart from this 
exception, the first order SW approximation is seen to approximate the exact results re- 
markably well even for small values of \a\ where quantum fluctuations evidently lead to a 
considerable reduction of M from its classical value. Notably, for negative values of a up 
to a ~ —0.1, the difference AM depends only weakly on the system size N. In this region, 
the second order correction M.2.sw{®) overestimates the difference between the first order 
SW approximation and the exact results. This is reminiscent of the SW expansion for the 
S = 1/2 square lattice HAF in which the third order correction is found to reduce the 
discrepancy between the exact results and the SW approximation by about 50% [[26J. For 



a — > ad, M.2.SW increases faster than AM. It is tempting to attribute this discrepancy to 
the unphysical divergence of M2.swi.Ci) at a c i. A selfconsistent SW treatment fT]| which 
eliminates this divergence is certainly called for. 



3.1b The collinear order 

As in the case of the 120° order we use the structure function Sjv(Q 2 ) at the Bragg 
vector Q2 of the collinear structure as a measure for the order parameter M . If the collinear 
order of the ground states of the finite cells were perfect the same steps that led to the 
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relations (FO) would yield 



WO = / -jf)- (3-5) 

The factor / counts the number of inequivalent Q2 vectors in the BZ. Hence, we get / = 3 
for the Cq v symmetric N = 12 cell, and / = 1 for the N = 18, 24, 30 cells which are not 
Cq v symmetric. 

As in the previous subsection we assume that the relation ( |3.5| ) still holds approximately 
when the ground state is not perfectly ordered. Figs. 7a,b,d and f show plots of Af 2 (a, N) 
for a>l/8 obtained by employing the numerical values of Sjv( Q2) in ( P-5|) . The maxima 
of J\f(a, 18) and AA(a,24) are seen to be only slightly reduced from the saturation value 
M = 1/2. This lends support to our assumption about the validity of ( |3.5| ). The stronger 
reduction of A/" 2 (a, 12) and the small discontinuity in A/" 2 (a,12) at a = 0.225 can be 
attributed to certain pecularities of the iV = 12 cell which will be discussed elsewhere |27] . 



In the transition region a ~ a c i, the behaviour of N{ot, N) is seen to differ qualitatively 
between the different system sizes. While J\f(a, 12) and 7V(a,24) grow smoothly with 
a, Af(a, 18) increases discontinuously at a qu (18) = 0.23. This discontinuity is due to a 
crossover in the ground state of the iV = 18 system whose translational symmetry changes 
from Q = to Q = Q 2 as a increases through a = a qu (18). In the classical picture 
( |1.2| ), the wavevector Q2 characterises the periodicity of the collinear state so that for 
a > a qu (18) the ground state of the iV = 18 system, IQ2), has the same periodicity 
as the classical state. For the N = 30 system we observe the same type of crossover 
between ground states at a qu (30) = 0.2. Because of limitations in computer memory 
we have not been able to calculate A/"(30) for a > a qu (30) where we would have had to 
work with a translationally noninvariant ground state |Q 2 ) - O ne can be sure, however, 
that the crossover between ground states will also result in a discontinuity in A/"(30) at 
ot = a qu (30). The difference between the iV = 12 and iV = 24 systems, for which the 
ground state remains at Q = for all values of a, and the iV = 18 and iV = 30 systems 
is reminiscent of Marshall's rule concerning the ground state of the square lattice HAF 
0: the ground state is translationally invariant, if the system size N is a multiple of 4; 
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otherwise, when N is even but not a multiple of 4, it has the periodicity of the classical 
state. 

As has been discussed in section 2 the SW theory suggests that the transition between 
the 120° structure and the collinear structure is of first order. The discontinuities in the 
order parameter M of the N — 18 and N = 30 systems directly confirm this prediction. 
For the system sizes iV = 12 and N = 24, the plots of Af and A4 are seen to overlap 
in a finite interval around a = a c i, Figs. 7a,d. While the overlap region decreases when 
the system size is doubled from iV = 12 to N = 24, the slopes of Af and A4 are seen to 
increase drastically at the transition points. These features agree with the behaviour of 
the order parameters of finite systems that exhibit a first order transition in the bulk limit. 
In summary, our numerical results support the view of a first order transition between the 
120° structure and the collinear structure of the HAFT. 




3.2 Helicity 

We wish to compare the SW results Xswi-N) and the exact values of the helicity. Numer- 
ically we calculate 

), (3-6) 

where Xa has been defined in ( |2.11| ). The identity ( |3.6| ) holds because the ground states 
are invariant under rotations in spin space. To find the connection between Xsw an d 
((Sa Xa) 2 ) we again consider the limit a — > — oo. Since the average ((J2a Xa) 2 ) contains 
four spin correlations an exact analytical evaluation is not possible in the limit. By a mean 
field approximation we find 

2. 




g S * + ^ + 0(N- 2 ), (3.7) 



a— »— oo 



mf 
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with a m f — ^ for S = ~. If Xi.sw/( a J N), ( |2.14j ), is evaluated for finite N one finds that 



the S 1 corrections vanish in the limit a — > — oo so that from ( p.14 ) and (|3] 




li ?M [xWKiV)] 2 = li ?co (f^i:^) ) - a -f + 0(N-' ) = 9S\ (3.8) 

Instead of working with the approximate finite size correction a m f we fit our exact results 
in the limit a ^ -oo by a second order polynomial in iV -1 , 

The fit is done separately for even N (N = 18, 24) and odd N (N = 21, 27). We find 

a>even = 7.63, b even = 20.02; a 0(M = 7.81, 6 odd = 5.75 . (3.10) 

Obviously, a m f = 7.75 is a good approximation for the leading finite size correction. That 
the finite size corrections of the helicity are large in comparison with the corresponding 
corrections of the sublattice magnetisation ( |3.3| ) is not surprising since each of the plaquet- 
tes in {(J2a Xa) 2 )o overlaps with six neighbouring plaquettes. In Figs. 9a-d we show plots 
of the numerical results 

^^((ip.))^-^ (3.11) 

and of the SW approximations x\ sw (a,N). As is seen in Fig. 10, where we plot the 
difference Ax(iV) = x(iV) — Xi.sw(N) between the first order SW results and the ex- 
act values for the helicity, the agreement is not as good as in the case of the sublattice 
magnetisation. The difference Ax(iV) deviates appreciably from the second order SW ap- 
proximation shown in Fig. 10 even for the large Cq v symmetric cells iV = 21 and iV = 27. 
Nevertheless, the enhancement of the helicity over the classical value, Figs. 9a-d, in the 
region — oo < a< — 0.2, which must be attributed to quantum fluctuations, is exhibited 
by both x{N) and Xi.Sw{N) in a similar fashion. Furthermore, the decrease of x(iV) and 
Xi.sw{N) to zero as a approaches a c i which is concomitant to the decrease of the sublattice 
magnetisation M. indicates clearly that the helicity is just another measure for the 120° 
long range spin order of the HAFT and has no independent significance as in the case of 
a spin nematic. 
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3.3 Chirality 

As we have discussed in the introduction, a RVB state which breaks time reversal and parity 
invariance has been proposed as the ground state of the HAFT ||18|| . Chiral symmetry 
breaking, i.e. a nonvanishing expectation value of the operator S$(Sj x Sjt) defined on an 
elementary triangular plaquette with sites i,j,k, is known to be the signature of such a 



state pqj . According to Baskaran |17j chiral symmetry breaking should only occur when 
the nnn coupling a exceeds a critical value. 
In Fig. 11 we plot the expectation value 

P ± = ±(C+C ± ), (3.12) 



where C± are the uniform and staggered chiralities 

2 

73 



C + = -7= E (Si(Si x S fe ) + S d (Si x S fc )), (3- 13a) 

£7 



C_ = v^E (S i (S i x Sfc) - S^S, x Sfc)). (3. 13b) 

zv 

Here, k and I denote the sites of a four-spin plaquette as depicted in Fig. 12, and the 
sums extend over the set of such plaquettes which covers the lattice. The normalisation in 
( j3. 134 b) is chosen such that the maximum modulus of the eigenvalues is unity in both 
cases. Contrary to Baskaran's conjecture and in agreement with the picture of the ground 
state as it has emerged so far no significant increase with a of either P + or P_ is observed. 



4 Summary and Conclusion 

In this work two important tools for investigating the ground state structure of quan- 
tum antiferromagnets, spin-wave theory and numerical diagonalisation of small periodic 
systems, have been applied to the HAFT with nn and nnn interactions. Our spin wave 
calculations extend previous work || on the 120° structure of the model in various ways. 
The extension of the S W-expansion of the sublattice magnetisation M. to order S~ 2 ap- 
pears to corroborate the previous conjecture |H| that the transition between the 120° order 
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and the collinear order of the HAFT is of first order. With the aim of comparing the SW 
results with the exact numerical results we develop finite size approximants for the former. 
By the example of the first order S W approximation for the sublattice magnetisation we 
demonstrate that the subleading 1/N size corrections must not be neglected for system 
sizes iV<10 2 . Moreover, irregular corrections which reflect the dependence of the finite size 
results on the shapes of the systems are seen to be relevant for smaller sizes iV<36. These 
results of the SW approximation indicate that extrapolations from system size iV<36 which 
exclusively rely on the leading 1/x/iV size dependence are unwarranted. Therefore, in the 
discussion of our exact numerical results we do not attempt to extrapolate to the bulk 
limit. Rather, we compare directly the finite size data obtained by the SW approximation 
with the corresponding exact data. 

In our numerical investigation we calculate quantities which are suited to characterise 
the 120° order and the collinear order of the ground state of the HAFT, namely, the order 
parameter M. and the helicity x y °f the 120° structure and the order parameter M of the 
collinear structure. To check whether there is an indication for unconventional order we 
also calculate the chirality P. Incommensurate structures which evolve in the classical 
picture for a > 1 are outside the scope of the present study since we work with finite 
periodic systems. 

For strong ferromagnetic nnn coupling, a <C — 1, one expects to find aggreement be- 
tween the exact results A4(a, N) and the first order SW approximation M.\.sw{.°<-> N) on 
physical grounds. Suprisingly, however, we find that M.i.sw( a i N) approximates Ji4(a, N) 
very well also for small values of \a\ where A4(a, N) is substantially reduced from its clas- 
sical value by quantum effects. This suggests that the SW approximation remains valid 
even for small positive values of a where the frustration already present in the triangular 
antiferromagnet is further enhanced by the nnn coupling. It appears that frustration is 
not as detrimental to long range order in spin-1/2 quantum antiferromagnets as has been 
thought [0, 13, 29]. The comparison of the first order SW approximation for the helicity 



X with the exact values supports this conclusion: \ is tied to the sublattice magnetisation 
just as it should be for an ordered 120° structure. 
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The behaviour of the order parameter Af(a, N) of the collinear structure is found to 
differ between the N = 12 and iV = 24 systems on the one hand and the N = 18 and 
N = 30 systems on the other. While Af(a, 12) and J\f(a, 24) increase smoothly as a in- 
creases through the classical transition point between the 120° order and the collinear order, 
Af(a, 18) and J\f(a, 30) increase discontinuously at a qu (18) = 0.23 and a qu (3Q) = 0.2 where 
A4(a, 18) and A4(a, 30) are still finite. These discontinuities are related to a crossover from 
a translationally invariant ground state for a < a qu to a ground state that has the peri- 
odicity of the collinear structure for a > a qu . Thus, for the system sizes N = 18 and 30 
one recovers the classical picture of a first order transition between the 120° order and 
the collinear order. In the cases N = 12 and 24 the ground states remain translation- 
ally invariant for all values of a. For these system sizes, the dependence of A4(a, N) and 
Af(a, N) on a and iV agrees with the finite size behaviour of the order parameters of a 
system that undergoes a first order transition between the ordered states characterised by 
Ai and Af in the bulk limit. In summary, our exact results support the conjecture [|l!| of 
a first order transition between 120° order and collinear order. In agreement with this we 
see no significant enhancement of the chirality in the transition region. 
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Appendix A: Second order spin wave approximation 



A straightforward evaluation of the expectation values in ( |2.8| ) yields the following expres- 



sions: 



where 



(H 2 )o 
N 



32S 2 



J> F(b) + -^--aG(b), 



1QS 2 



B 
3S 



F(b) = i [ 2-2J 1 + -^(b-2a(k))(A 1 (k,B) + A 2 (k,B))j 

+Bl + AB 1 B 2 - 2B\ - {B 1 + 2B 2 )(A X - A 2 ), 
G{b) = (A 1 + A 2 -2-C 1 -C 2 ) 2 + (C 1 -C 2 )(C 1 -C 2 -A 1 + A 2 ) 



(A.l) 



and 



with i/(k, B) as defined in (O). Ai, Bi and Ci denote the sums 



where 



Bi 

a 



iv k 



N 



ET(k)A(k), 



i = 1,2 



Ai(k) = [A 2 (k)] 



-i 



1 + b - 2a(k) + 2 7 (k) 



\J 1 + 6 - 2a(k) - 7(k) 

with a(k)) and 7(k)) as defined in ( |2. 6a| ,b). 
Similarly, 



E 



(K}|#3/ 2 |K}> 



3 1 



E 



/ 2 (k, P ) 



{ £ } £({0}) - £({n k }) 16S 2 iV 2 £ Kk) + Kp) + Kk + P 



(A.2) 
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where 



/ 2 (k )P ) 



with 



0(k) 



^(k^p^k + p) [/3(k)(l - A 2 (p)A 2 (k + p)) 
+/3(p)(l-A 2 (k)A 2 (k + P )) 
-/3(k + p)(l-A 2 (k)A 2 (p))] 2 , 



• /7 \ • / 1 r \/3, . . ,1 , v^3 , s 
sm(Kj - sm(-fc,. - -g-fey) - sin^-fc, + —k y ) 



Miyake's result || is reproduced by setting a = in these expressions. 

To obtain the second order correction Xz.sw to the helicity we follow the same procedure 
as in the determination of M.2.SW- We replace the Zeeman term in (|2.1| ) by the term 



y 

A) 



which couples the helicity to a fictitious source field h so that 

1 d 



X 



y — — 



Ndh 



E Q (a,h) 



where E (a, h) is the /i-dependent ground state energy of the system. The second order 
term in the expansion of E (a, h) with respect to S^ 1 has the same structure as E (a, B), 
( p.8|) , i.e. it consists of two additive contributions that contain singular parts proportional 
to \fh which cancel in the sum. Thus, the limit 

EP(a,h)-EP(a,0) 



xt.sw 



lim 

N 



h 



can be performed numerically in exactly the same way as in (|2.10|) . 



Appendix B: M-i.sw m the limit —ol/N — > oo 



From (|2~6a| ,b) it follows that for k ^ 0, ±Qi 



a 



l«(k)i>37r 2 y 



(b.i; 
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and hence, in the limit \a\/N ^> 1, 

|a(k)|»l>| 7 (k)|, 
so that Ai t 2(k) can be expanded in this limit: 

^(k) = l T i±^ffl ± lf# + Q(JV/^. 



2a(k) 2a(k) 
Using this expansion in ( 2.17|) one finds the physically correct result 



Mi. sw (a,N) = S + 0(N 2 /a 2 ) 



By contrast, the expression 



M 1 .sw(<x,N) = S{l + 



4S 



2 4e e mv) 

i iV KeSZ iV KeSZ 



(B-2) 



(B-3) 



(B.4) 



(B.5) 



which derives from ( |2.7|) if one neglects only the infinite contributions of the zero modes, 
yields the unphysical result 

1 6 



M 1 . sw {a, N) = S { 1 + — ^ + 0(N 2 /a 2 ' 



(B.6) 



Obviously, the finite size correction —3 / (2N) which is accounted for in ( 2.171) is important 
for systems sizes N < 36 for which exact numerical results are available. 



Appendix C: Connection between the structure func- 
tion and the order parameter 



To evaluate SV(Qi) in the limit —a > 1 we start from the expression 



/ /N/3 N/3 \ v 

sw(Qi) = 3 ( E - E Us,) 



(C.l) 



which follows from (|3.1| ) by taking the equivalence of the three sublattices into account. 
Since the spins of one sublattice are ferromagnetic aligned for —a > 1 we have 



lim (SjSj) = S 2 for i,jeA, i ^ j. 



a— >— oo 



(C.2) 
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However 



lim (SSj) = S 2 c for zeA, jeB, (C.3) 



a— >— oo 



where c can differ from the value cos(27r/3) = —1/2, since the spins of different sublattices 
may fluctuate against each other, c is fixed by the condition that the ground state of 
systems of size N must be a singlet or a doublet depending on whether N is even or odd: 

N 

lim S tot (N) = lim (QTS;) 2 ) 

i 

N 2 ( 6 \ , iVeven, 

= v? 1 + 2c+ 7v = , (a4) 



4 ' 



Hence, 



, A^even, 
-- — — 4- — A^ndd 

2 Af 2iV 2 ' UUU - 



(C.5) 



Using ( |C.2| ), ( |C3| ) and ( |C5| ) we find from ( |C1| ) the expressions ( |3.3| ) of the main text for 

^(Qi). 
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Figure Captions 



Fig.l Classical ground state configuration in the regime a < 1/8. Bold lines: nn 
bonds; dashed bold line: nnn bond 

Fig. 2 Classical ground state configuration in the regime 1/8 < a < 1. The state 
is degenerate with respect to the angle 9. 

Fig. 3 The Brillouin zone of the triangular lattice. The arrows indicate the path 
which the Bragg vector Q of the ground state configuration takes when 
a increases: Q = Qi for — oo < a < 1/8; Q = Q 2 for 1/8 < a < 1; 
for 1 < a < oo Q moves continuously from Q 2 towards the corner of the 
v^x^ Brillouin zone. 

Fig. 4 a,b First and second order spin wave results: (a) the sublattice magnetisation 
M. as a function a; (b) the helicity x y as a function of a. 

Fig. 5 Size dependence of the first order spin wave approximation for M. at a = 0. 

: -M.i.sw( a — 0,iV); full line: leading size correction to M.\.sw{^i °°); 
dashed line: leading and subleading correction to AAi.swi®, oo) 

Fig. 6 a-f Finite systems studied in this work. For the N = 18, 24 and 30 cells which 
are not C 6v symmetric, the triangular lattice has been distorted such that 
the outer contours of the cells take the shape of regular hexagons. The 
degree of distortion illustrates the deviation from the Cq v symmetry. 
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Fig. 7 a-f The square of the order parameters A4 and J\f of the 120° structure and of 
the collinear structure as a function of a for various system sizes. M.\.sw 
is obtained from eq. fl2.17|) . The insets in (a) and (d) show an enlargement 



of the transition region. 

Fig. 8 The difference AM{a, N) = M{a, N)—Mi.sw{®, N) and the second order 
spin wave correction M.2.sw(p) as functions of a. 

Fig. 9 a-d The square of the helicity x y °f the 120° structure as a function of a for 
various system sizes. Xi.sw is obtained from eq. ( p.21| ). 

Fig. 10 The difference Ax(a,iV) = x( a ,N) — Xi.sw( a , N) and the second order 
spin wave correction X2.swi.Ci) as functions of a. 

Fig. 11 The chiralities P±(a,N) as functions of a. 

Fig. 12 A plaquette of four spins on which the uniform and staggered chiral opera- 
tors are defined. 
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